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Abstract 

We present measurements of CP-violation parameters in 6 ^ 57 transitions based on a sample of 
386 X 10^ BE pairs collected at the T(4S') resonance with the Belle detector at the KEKB energy- 
asymmetric e+e^ collider. One neutral B meson is fully reconstructed in the B^ K^tt^^ decay 
channel irrespective to the K^ir^ intermediate state. The flavor of the accompanying B meson is 
identified from its decay products. CP-violation parameters are obtained from the asymmetries in 
the distributions of the proper-time intervals between the two B decays. 

We obtain the following results for the K^gir*^ invariant mass covering the full range up to 
l.SGeV/c^: 

5^o^o^ = -K0.08±0.41(stat) ±0.10(syst), 
^^0^0^ = -K0.12±0.27(stat) ±0.10(syst). 

PACS numbers: ll.30.Er, 12.15.Hh, 13.25.Hw 
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INTRODUCTION 



In the Standard Model (SM), CP violation arises from an irreducible phase, the 
Kobayashi-Maskawa (KM) phase [1], in the weak-interaction quark-mixing matrix. The 
phenomena of time-dependent CP violation in decays through radiative penguin processes 
such as 6 S7 are sensitive to physics beyond the SM. Within the SM, the photon emitted 
from a {B^) meson is dominantly right-handed (left-handed). Therefore the polarization 
of the photon carries information on the original 6-flavor and the decay is, thus, almost 
flavor-specific. As a result, the SM predicts a small asymmetry [2, 3] and any significant 
deviation from this expectation would be a manifestation of new physics. It was pointed 
out that in decays of the type B^ P^Q^'-j, where P° and represent any CP eigen- 
state spin-0 neutral particles (e.g. P° = and Q° = 7r°), many new physics effects on 
the mixing-induced CP violation do not depend on the resonant structure of the P^Q^ 
system [4]. 

At the KEKB energy-asymmetric e+e~ (3.5 on 8.0 GeV) collider [5], the T(4S') is pro- 
duced with a Lorentz boost of j3'^ = 0.425 along the z axis, which is defined as the direction 
antiparallel to the e"*" beam direction. In the decay chain T(4S') B^B^ — > /sig/tag, where 
one of the B mesons decays at time tgig to a final state /gig, which is our signal mode, and 
the other decays at time ttag to a final state /tag that distinguishes between B^ and B^, the 
decay rate has a time dependence given by 

p-i^'i/^B" r r ■^^ 

V{At) = — <^ 1 + q Ssm{AmdAt) + Acos{AmdAt) }. (1) 

Atbo L J J 

Here S and A are CP- violation parameters, r^o is the P" lifetime, Am^ is the mass difference 
between the two P° mass eigenstates. At is the time difference 4ig — ^tag, and the 6-fiavor 
charge q — +1 (—1) when the tagging B meson is a P° (P"). Since the P° and P° mesons 
are approximately at rest in the T(4S') center-of-mass system (c.m.s.). At can be determined 
from the displacement in z between the /gig and /tag decay vertices: At ~ (-^sig^-^^tag) / if^ic) = 
Az/iP^c). 

For B^ — > KgTr'^j, the Kg vertex is displaced from the B vertex and often lies outside 
of the silicon vertex detector (SVD). When the Kg vertex can be reconstructed inside the 
SVD, the time-dependent CP asymmetry can be measured. Measurements of such CP 
asymmetries were previously reported by BaBar [6] and Belle [7] in the P° K*^{^ 
KgTr'^)^ decay: 

Sk*o^ = 0.25 ±0.63 ±0.14 (BaBar) 
Sk^o^ = -0. 7912 ±0.10 (Belle). 

Belle also measured these asymmetries with an extended M^o^o mass region [7] {Mj^o^o < 
1.8GeV/c2): 

SKO^o^^-0.58t',:tl±0.n (Belle). 

In this analysis, we update the CP measurements for P'' KgTT^'-f in the mass region 
^K^wO < 1.8GeV/c^ with an additional dataset of 111 x 10^ PP pairs. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of an SVD, 
a 50-layer central drift chamber (CDC), an array of aerogel threshold Cerenkov counters 
(ACC), a barrel-hke arrangement of time-of-fiight scintillation counters (TOF), and an 



4 



electromagnetic calorimeter comprised of CsI(TI) crystals (ECL) located inside a super- 
conducting solenoid coil that provides a 1.5 T magnetic field. An iron flux-return located 
outside of the coil is instrumented to detect mesons and to identify muons (KLM). The 
detector is described in detail elsewhere [8] . Two inner detector configurations were used. A 
2.0 cm beampipe and a 3-layer silicon vertex detector (SVDl) was used for the first sample 
of 152 X 10^ BB pairs, while a 1.5 cm beampipe, a 4-layer silicon detector (SVD2) and a 
small-cell inner drift chamber were used to record the remaining 234 x 10^ BB pairs [9]. 

EVENT SELECTION, FLAVOR TAGGING AND VERTEX RECONSTRUCTION 
Event Selection for K^-k^^ 

For high energy prompt photons, we select an isolated cluster in the ECL that has no 
corresponding charged track, and has the largest energy in the c.m.s. We require the shower 
shape to be consistent with that of a photon. In order to reduce the background from 7r° and 
r] mesons, we exclude photons compatible with 7r° — > 77 or — > 77 decays; we reject photon 
pairs that satisfy Ct,o > 0.18 or £^ > 0.18, where £7rO(r?) is a tt" (77) likelihood described in 
detail elsewhere [10]. The polar angle of the photon direction in the laboratory frame is 
restricted to the barrel region of the ECL (33° < 9^ < 128°), but is extended to the end-cap 
regions (17° < 9^ < 150°) for the second data sample due to the reduced material in front 
of the ECL. 

Neutral kaons {Kg) are reconstructed from two oppositely charged pions that have an in- 
variant mass within ±6MeV/c^ {2a) of the Kg nominal mass. The 7r"'"7r~ vertex is required 
to be displaced from the interaction point (IP) in the direction of the pion pair momen- 
tum [11]. Neutral pions {n^) are formed from two photons with the invariant mass within 
±16MeV/c^ (30") of the 7r° mass. The photon momenta are then recalculated with a vr'' 
mass constraint and we require the momentum of 7r° candidates in the c.m.s. to be greater 
than 0.3GeV/c. The K^n'^ invariant mass, M^o^o, is required to be less than 1.8GeV/c^. 

B^ mesons are reconstructed by combining Kg, 7r° and 7 candidates. We form two 
kinematic variables: the energy difference AE = — E^-^^- and the beam-energy 

constrained mass Mbc = ^(-^be^m)^ ~ (Pb°^'^')^' where i?bcam' the beam energy, i?^™-^- 
and p^'"^ *^- are the energy and the momentum of the candidate in the c.m.s. Candidates are 
accepted if they have Mbc > 5.2 GeV/c^ and -0.5 GeV < AE < 0.5 GeV. 

We reconstruct S+ K^n'^^ candidates in a similar way as the B° — > K^n^^ decay 
in order to reduce the cross-feed background from B~^ KgH^j in B^ K^n^^. The 
B~^ Kgii^^ events are also used for various crosschecks. For a 7r+ candidate, we require 
that the track originates from the IP and that the transverse momentum is greater than 
0.1 GeV/c. We also require that the tt"*" candidate cannot be identified as any other particle 
species (ii'+,p+, e"*", and ^~^). 

Candidate B~^ KgH^'-f and B^ Kgii^^ decays are selected simultaneously; we allow 
only one candidate for each event. The best candidate selection is based on the event 
likelihood ratio T^-g/b that is obtained from a Fisher discriminant J-" [12], which uses the 
extended modified Fox- Wolfram moments [13] as discriminating variables. We select the 
candidate that has the largest TZs/h- The signal region is defined as —0.2 GeV < AE < 
0.1 GeV and 5.27GcV/c2 < Afbc < 5.29GeV/c2. 

We use events outside the signal region as well as large Monte Carlo (MC) samples 
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to study the background components. The dominant background is from continuum hght 
quark pair production (e+e~ ^ qq with q = u,d,s,c), which we refer to as qq hereafter. 
In order to reduce the qq background contribution, we form another event hkehhood ratio 
^?/b combining TZs/h with cos9h and cos^b, where is the polar angle of the B meson 
candidate momentum in the laboratory frame, and 6h is the angle between the B candidate 
momentum and the daughter Kg momentum in the rest frame of the K^ti system. Since 
the relative background contribution will be smaller in the region of the K*, we introduce 
two KgTr'^ invariant mass regions: MRl, defined as 0.8GeV/c^ < M^o^o < l.OGeV/c^, 
and MR2 which is defined as M^o^o < l.SGeV/c^ after excluding MRl. The specific TZfJ^ 
selection criteria applied depend on both the mass region and flavor tagging information. 
After applying all other selection criteria described so far, 77% of the qq background is 
rejected while 87% of the K*'^j signal is retained in MRl; in MR2, 87% of the qq is rejected 
while 68% of the K2^'y signal is retained. Background contributions from B decays, which 
are considerably smaller than qq, are dominated by cross-feed from other radiative B decays 
including S+ — > K^n'^^. 

Flavor Tagging 

The 6-fiavor of the accompanying B meson is identified from inclusive properties of par- 
ticles that are not associated with the reconstructed signal decay. The algorithm for flavor 
tagging is described in detail elsewhere [14]. We use two parameters, q defined in Eq. (1) and 
r, to represent the tagging information. The parameter r is an event-by-event flavor-tagging 
dilution factor that ranges from to 1; r = when there is no flavor discrimination and 
r — 1 implies unambiguous flavor assignment. It is determined by using MC data and is 
only used to sort data into six r intervals. The wrong tag fraction w and the difference 
Aw in w between the S° and S° decays are determined for each of the six r intervals from 
data [11]. 

Vertex Reconstruction 

The vertex position of the signal-side decay is reconstructed from the Kg trajectory with 
a constraint on the IP; the IP profile {ax ^ 100 //m, ay ~ 5 /im, cr^ ~ 3 mm) is convolved 
with the finite B flight length in the plane perpendicular to the z axis. Both pions from 
the Kg decay are required to have enough hits in the SVD in order to reconstruct the Kg 
trajectory with high resolution: at least one layer with hits on both sides and at least one 
additional hit in the z side of the other layers for SVDl, and at least two layers with hits on 
both sides for SVD2. The reconstruction efficiency depends not only on the Kg momentum 
but also on the SVD geometry. The efficiency with SVD2 (51%) is significantly higher than 
with SVDl (40%) because of the larger detection volume. The other (tag-side) B vertex 
determination is the same as that for the B^ — > (pKg analysis [11]. 

SIGNAL YIELD EXTRACTION 

Figure 1 shows the Mbc (AE) distribution for the reconstructed Kgir^^ candidates within 
the AE (Mbc) signal region after flavor tagging and vertex reconstruction. The signal yield 



6 




M^^ (GeV/c^) AE(GeV) 

FIG. 1: (a) Mbc distributions within the AE signal region and (b) AE distributions within the Mbc 
signal (MRl and MR2 combined). Solid curves show the fit to signal plus background distributions. 
Lower dashed curves show the background contributions from qq. Upper dashed histogram show 
the sum of background contributions from qq and B decays. 

is determined from an unbinned two-dimensional maximum-likelihood fit to the AE-My,c 
distribution. The fit region is chosen as -0.4 GeV < AE < 0.5 GeV and 5.2 GeV/c^ < Mbc 
to avoid other BB background events that populate the low-Ai? high-M^^o^o region. The 
signal distribution is represented by a PDF obtained from an MC simulation of B^ K*^'y 
and B^ -ft'j'^T i^d^t accounts for a small correlation between Mbc and AE. The background 
from B decays are also modeled with an MC simulation. For the qq background, we use 
the ARGUS parameterization [15] for Mbc and a second-order polynomial for AE. The 
normalizations of the signal and background distributions and the qq background shape 
are the five free parameters in the fit. We observe a total of 260 candidates in the signal 
box in MRl, which decreases to 116 after flavor tagging and B vertex reconstruction, and 
obtain 70 ± 11 signal events from the fit; the average signal purity over the six r intervals is 
63 ± 11%. In MR2, corresponding numbers are 236, 120, 45 ± 11, and 41 ± 11%. 

CP ASYMMETRY MEASUREMENTS 

We determine S and A from an unbinned maximum-likelihood fit to the observed At 
distribution. The probability density function (PDF) expected for the signal distribution, 
Vsig{At;S,A,q,w, Aw), is given by the time dependent decay rate [Eq. (1)] modified to 
incorporate the effect of incorrect flavor assignment. The distribution is convolved with 
the proper-time interval resolution function -Rsig, which takes into account the flnite vertex 
resolution. The parametrization of i?sig is the same as that used for the B^ ^s'^^ 
decay [11]. i?sig is flrst derived from flavor-speciflc B decays [16] and modifled by additional 
parameters that rescale vertex errors to account for the fact that there is no track directly 
originating from the B meson decay point. 
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For each event, the following likehhood function is evaluated: 



Pi -(1 - /ol) / /sigPsig(At')i?sig(At, - M') 



J —oo 



+ (1 - /sig)nkg(At')i?bkg(Ati - At') d(At') 



(2) 



+ folPol{Ati), 



where Poi is a Gaussian function that represents a small outlier component with fraction 
/ol [17]. The signal probability /sig is calculated on an event-by-event basis from the func- 
tion which we obtained as the result of the two-dimensional AE-M]^^ fit for the signal yield 
extraction. A PDF for background events, Vhkg, is modeled as a sum of exponential and 
prompt components, and is convolved with a Gaussian which represents the resolution func- 
tion i?bkg for the background. All parameters in Vhkg and Rhkg are determined by a fit to the 
At distribution of a background-enhanced control sample, i.e. events outside of the AE'-M^c 
signal region. We fix r^o and Am^ at their world-average values [18]. 

The only free parameters in the final fit are Sj^a^o^ and w4j^o^o^, which are determined 
by maximizing the likelihood function L = Yli Pi{A.ti; S , A) where the product is over all 
events. We obtain 



We define the raw asymmetry in each At bin by [Ng=^i — Ng=^i) / {Ng=^i + Nq=„i), 
where A''g=_|_i(_i) is the number of observed candidates with q = +1(— 1). Figure 2 shows 
the raw asymmetries for the K^Ti'^'y events. Note that these are simple projections onto the 
At axis, and do not reflect other event-by-event information (such as the signal fraction, 
the wrong tag fraction and the vertex resolution), which is in fact used in the unbinned 
maximum-likelihood fit for S and A. 

Systematic Error 

Primary sources of the systematic error are (1) uncertainties in the resolution function 
(±0.06 for S^Ot^Oj and ±0.03 for ^j^o^o^), (2) uncertainties in the vertex reconstruction 
(±0.03 for iSj^o^o^ and ±0.04 for ^^o^o^) and (3) uncertainties in the background fraction 
(±0.07 for tS^o^o,, and ±0.03 for ^^o^o^). Effects of tag-side interference [20] contribute 
±0.07 for Ax^^tt"^-- Also included are effects from uncertainties in the wrong tag fraction and 
physics parameters {Ama, tbo and Ak'Oj). Fitting a large sample of MC events revealed no 
bias in the fit procedure. The statistical errors from the MC fit are assigned as systematic 
errors. The total systematic error is obtained by adding these contributions in quadrature. 

Crosschecks 

Various crosschecks of the measurement are performed. We apply the same fit proce- 
dure to the 5° — > J/ipKg sample without using J/ip daughter tracks for the vertex recon- 



struction [21]. We obtain Sj/^k^ = ±0-73 ± 0.08(stat) and Aj/^k^ = ±0-01 ± 0.04(stat), 



±0.08 ± 0.41(stat) ± O.lO(syst), 
±0.12 ± 0.27(stat) ± O.lO(syst). 



8 




-2.5 2.5 
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-2.5 2.5 
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FIG. 2: Asymmetry in each At bin for S° K^tt'^j with < r < 0.5(left) and 0.5 < r < l.O(right). 
Sohd curves show the results of the unbinned maximum-hkehhood fits. 



which are in good agreement with the world- average values [19]. We perform a fit to 
KgTc'^'-f , which is a good control sample of the Kgir^'j decay, without us- 

ing the primary vr"*" for the vertex reconstruction. The result is consistent with no CP 
asymmetry, as expected. Lifetime measurements are also performed for these modes, and 
values consistent with the world-average values are obtained. Ensemble tests are car- 
ried out with MC pseudo-experiments using S and A obtained by the fit as the input 
parameters. We find that the statistical errors obtained in our measurements are all 
within the expectations from the ensemble tests. Fits to the two M^o^o regions yield 
S = -fO.Ol ± 0.52(stat) ± O.ll(syst) and A = +0.11 ± 0.33(stat) ± 0.09(syst) for MRl, and 
S = +0.20 ± 0.66(stat) and A = +0.14 ± 0.46(stat) for MR2. The results are consistent 
with those from the full Mj^o^o sample. 



SUMMARY 

We have performed a measurement of the time-dependent CP asymmetry in the decay 
B^ — > Kgii^'-f with KgTT^ invariant mass up to 1.8 GeV/c^, based on a sample of 386 x 10^ BB 
pairs. We obtain CP-violation parameters 5^0,^0^ = +0.08 ± 0.41(stat) ± O.lO(syst) and 
•^K^-n-o-y = +0.12 ± 0.27(stat) ±0.10(syst). We do not find any significant CP asymmetry, 
and therefore no indication of new physics from right handed currents, with the present 
statistics. 
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